Charge Fluctuations in a Quantum Dot with a Dissipative Environment 
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We consider a multiple tunneling process into a quantum dot capacitively coupled to a dissipa- 
tive environment. The problem is mapped onto an anisotropic Kondo model in its Coulomb gas 
representation. The tunneling barrier resistance and the dissipative resistance of the environment 
correspond to the transverse and the longitudinal Kondo couplings respectively. We thus identify 
a line in the parameter space of the problem which corresponds to a zero-temperature Berezinskii- 
Kosterlitz-Thouless like phase transition. The physics of coupling to the environment is elucidated 
and experimental consequences of the predicted transition are discussed. 
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I. INTRODUCTION 

The problem of the Coulomb blockade in zero dimen- 
sional quantum dots has been introduced theoretically 
and studied both theoretically and experimentally over 
the past decade [[l]j|] . First studies employed an orthodox 
model within which the effect of charging has been ac- 
counted for classically. More recent studies have included 
various quantum mechanical aspects of the Coulomb in- 
teraction on various levels of rigor. Thus the Altshuler 
Aronov zero bias anomaly || has been related to the 
Coulomb blockade [jj. In a slightly different language 
it was also realized that the coupling to a dynamically 
polarizable environment may modify the system's behav- 
ior qualitatively ■ The instantaneous tunneling of an 
electron into the dot leads (through the Coulomb inter- 
action) to a shake-up excitation of the low-energy modes 
of the environment. This process, which is very similar in 
nature to the X-ray edge singularity (see Refs. J||| an d 
references therein) , results in a power law current -voltage 
characteristics of a single tunneling barrier connected to 
a dissipative circuit. 

Another remarkable sequence of developments is as- 
sociated with the study of the quantum fluctuation of 
the charge on the dot, or multiple tunneling processes 
[p~0| — |l6|| . The importance of such processes in the vicinity 
of transmission resonance has been first realized by Glaz- 
man and Matveev [jy]. Later on Matveev |Q mapped 
the problem onto the anisotropic Kondo model |]l7| , |l8| . 
This mapping is based on a projection of the Hamilto- 
nian onto two charge states of the dot. The two degen- 
erate charge states of the dot (which corresponds, say, to 
N and N + 1 electrons respectively), are mapped onto 
two spin states, "up" and "down" , with the coupling to 
the leads playing the role of a coupling to a Fermi sea. 
A small deviation from resonance conditions plays the 
role of a constant magnetic field in the Kondo problem, 
leading to Zceman splitting. Extensions of this problem 
to other variants have been carried out employing equiv- 
alent methods [fill. 



In the present paper we consider the effect of a dissipa- 
tive environment on the quantum fluctuations of charge 
of the dot. This relatively straightforward combination 
of previously studied effects leads to a qualitatively novel 
physical situation, which, to the best of our knowledge, 
has not been discussed before in the present context. The 
dissipative resistance of the environment, r, capacitively 
coupled to the dot, gives rise to a new parameter, im- 
portant for the physical description of the system. In the 
language of the anisotropic Kondo model Jlj],[l8| this new 
parameter is precisely equivalent to a (ferromagnetic) z- 
axis coupling, J z . Let us recall that the resistance of the 
tunneling barrier, r t , is mapped |l2| ] onto the (inverse) x- 
y coupling, J±. As a result one obtains a separatrix line 
in the parameter space (r t -r) of the system. Crossing 
this line by varying the values of the respective resis- 
tances leads to a zero-temperature phase transition. In 
the "ferromagnetic" phase the system is overdamped and 
its behavior is akin to that of a single tunneling barrier 
connected to a dissipative circuit |^,^] . Namely, the peak 
conductance decreases as a certain power of temperature. 
By contrast, in the "anti-ferromagnetic" phase, the peak 
hight undergoes re-entrance behavior and eventually ap- 
proaches a perfect conductance limit |Q (for a single 
channel tunneling contact). 

Technically we proceed by mapping the problem of 
multiple tunneling into the dot, onto that of a classical 
neutral gas in one dimension with long range interac- 
tions. In the past, such a mapping has been employed 
in the context of the anisotropic Kondo problem [p7|Jl8[ 
and dissipative quantum tunneling |pOf . In the latter 
context the analogous phase transition corresponding to 
localization-delocalization in double quantum well was 
exploited extensively. The quantum dot realization of 
this transition, proposed here, may be the most suitable 
for experimental realization. A related, though quali- 
tatively different gas, which also can be realized in the 
present context, has been employed in the study of a 
one-dimensional Luttinger model with an impurity pi) . 
Here we restrict ourselves to a study of a quantum dot 



1 



connected to a lead through a single tunnel barrier, occa- 
sionally referred to in the literature as a "single electron 
box" . We discuss the mapping onto a classical gas and 
the subsequent renormalization group (RG) procedure 
on the level of the partition function of such a "box" . A 
quantity more directly related to the experiment is the 
conductance through a dot with the two barriers jl5],|l6| ■ 
We believe, though, that our calculations reflect the ex- 
perimental situation with rt being the resistance of the 
weakest of the two barriers. We also require that our dot 
is large enough, in a sense that one may disregard the 
discreteness of its single particle spectrum. 

The outline of this paper is the following. In Section |l| 
we define our problem in terms of a model Hamiltonian 
and an equivalent electrical circuit, and discuss qualita- 
tively the effect of the environment. We also discuss a 
complete solution of the problem in lowest order in the 
tunneling amplitude. Mapping the quantum fluctuations 
of the charge onto a classical Coulomb gas is presented 
in Section III. The fully and partially screened cases are 



solved using a renormalization group procedure. In Sec- 
tion |y| we discuss the physics of our results and briefly 
comment on possible generalizations. Technical details 
of the derivation are presented in two Appendices. 



the bare (high temperature) tunneling resistance, Rt 
(27rh/e 2 )r t , where 



(rt)- 1 = (2nf\W\ 2 v 



[o] v m 



(4) 



where z/M and lM are the (bare) density of states (DOS) 
of the dot and the lead respectively. Finally, we shall 
assume that the dot is embedded in an electromagnetic 
environment, felt by the dot as an external noise. This 
noisy background is described as an effective time depen- 
dent term in the Hamiltonian 



Tin 



= »?(*)£« 



(5) 



Below we specify the noise, r)(t), in terms of the environ- 
ment's characteristics. Equivalently, we could say that 
the dot is connected to a bath of harmonic oscillators 
which are integrated out. The spectral density of r)(t) re- 
flects the nature of the environment |2(J . In this sense the 
averaging over noise should be understood as annealed 
and not as quenched. The total effective Hamiltonian of 
the system is thus 



H = TLdo 



H 



lead 



K 



coupling 



(6) 



II. MODEL HAMILTONIAN AND AN 
EQUIVALENT CIRCUIT 

We consider a small capacitance dot, described by the 
Hamiltonian 



Ho — e a a<a a a 




(1) 



where a^(a a ) is a creation (annihilation) operator of an 
electron in an exact single particle state a (the latter is 
defined including disorder potential and spin); V^/2 = 
e 2 /2C is the charging energy, associated with the total 
charge on the dot (C is the dot's self-capacitance); Afo is 
the effective charge of the positive background. This may 
be regarded as the large wave-length (q=0) component 
of the interaction. Other, q =/= 0, components are ignored. 
The dot is assumed to be weakly coupled to a conducting 
lead, described by a gas of non-interacting electrons 



H 



lead 



k 



(2) 



where dk,d~£ are Fermi operators. The coupling to the 
lead is given by 



K 



coupling / 
fc , a 



(W ka d+a a + W£ a a+d k ) 



(3) 



Below we shall consider a single point contact and im- 
pose Wka — W p2|. The parameter W is related to 



The physical effect of the environment |q|| can be un- 
derstood in the following manner. Imagine that an elec- 
tron has been injected into the dot. The total effective 
charge, e((N) — A/o), interacts with charges in the en- 
vironment, leading to their redistribution (polarization). 
The polarized charge of the environment reduces the en- 
ergy cost of adding (or removing) an electron to (or from) 
the dot. There is a certain finite time scale characterizing 
this redistribution of the environmental charge, rendering 
the effective interaction in the dot non-instantaneous (re- 
tarded) [p3p|] . This time constant can be described as an 
RC time of an effective electric circuit. It might appear 
that the results obtained through the ensuing analysis 
are highly non-universal and depend on the particular 
choice of the model for the environment. We stress that 
the effect of the environment is quite general, the only 
model dependent feature being the concrete form of the 
screened zero-mode interaction, V(u>) (see below). It 
enters through the dependence of V(uj) on an effective 
impedance of the environment |24{| , and is due to the 
fact that details of the environment's polarization may 
depend on various objects located far from the dot. 

Fig. ^ shows a simplified equivalent circuit depicting 
an experiment of an electron injection into a quantum 
dot. The dot is assumed to be weakly coupled to the 
lead; C represents the self-capacitance of the dot. The 
dot is capacitively coupled (through Cg) to a dissipative 
environment whose impedance is Z(oj). The source U 
represents the equilibrium voltage noise of the environ- 
ment. 
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FIG. |IJ An equivalent circuit describing an experimental 
setup of an electron injection into a dot. For details see the 
text. 



In the lowest order in the tunneling amplitude, \W\ 2 , 
the conductance of the tunneling junction is fully deter- 
mined by the DOS of the dot. The latter may be obtained 
through the imaginary-time one particle Green function 
of the Hamiltonian TIq + Tl no i se . It is, of course, strongly 
modified by the effects of interaction and coupling to the 
environment. For completeness we briefly reproduce the 
calculations of Ref. in Appendix A, which lead to the 
following expression for the one particle Green function 



-S(r) 



Here the effective action, S, is given by 



S(t) = T 



(7) 



(8) 



m/0 



where T is the temperature and io m are the bosonic 
Matsubara frequencies, Lu m = 2i:mT. The form of the 
screened interaction, V(u> m ), can be easily deduced rely- 
ing on the equivalent circuit picture, Fig. [l]. It is shown 
in Appendix A that 



V(u m ) = V {0) - K(u m ) 



(9) 



where K is the correlator of the Gaussian noise felt 
by a dot. In the imaginary time domain K(t — 
T ') = { 7 l{ T ) T l{ T ')) i an d may be evaluated employing the 
fluctuation-dissipation theorem. Consulting Fig. |l| the 
correlator of the total noise in the equivalent circuit (in 
terms of real frequencies) is: (UU)^ — e 2 iu]Zt t{u) , 
where the total impedance of the equivalent circuit is 
Z to t — (iuC)^ 1 + (iujCg)^ 1 + Z(u) . The corresponding- 
voltage drop on the dot is ij — U^ujC)^ 1 /Z to t ■ Substi- 
tuting the above given expressions in Eq. (^|) and rewrit- 
ing it in a finite temperature form one obtains 



C/C c 



\zc 



C l + C/C G + \Lu m \ZC 



(10) 



In the high frequency limit the interaction is unscreened, 
whereas in the low frequency limit we obtain a (partially) 



screened interaction, where C is replaced by C + Cq- 
The scale of the crossover frequency is given by Vl^ 1 = 

iz^icn. 

Coming back to the action, Eq. (||) , the long time be- 
havior of S(t) is given by the small frequency limit of 
V(u m ), yielding 



S(t) 



1 2(C + C G ) 



for ZC G <r</3 



(11) 



This is the case of a partially screened interaction, C — > 
(C + Cg). One expects that in this case the addition of 
charge to the dot costs finite energy even at long times, 
which is why Eq. ([n]) is linear in r. This immediately 
implies the existence of a gap in the DOS, which is a 
manifestation of the Coulomb blockade. In the special 
case of a fully screened interaction (Cq — ► oo), the long 
time interaction vanishes, i.e. we expect V(u m — > 0) = 0. 
This is the case of a single tunneling contact in series with 
a classical impedance, Z, considered in Refs. [§,[§. F° r a 
purely Ohmic environment, Z = i?, the interaction is 



V{u m ) = 



(12) 



with fl = (RC)^ 1 . On short time scales (u> > f2) before 
the environment becomes fully polarized to screen out the 
extra charge injected into the dot, the charging energy is 
finite. Thus the injection of an electron can be considered 
as "tunneling under the barrier" in the time direction. 
The very same physics has been discussed in the context 
of tunneling into two dimensional systems with diffusive 
disorder |2^,^|. The related energy cost on short time 
scales suppresses free particle exchange between the dot 
and the particle reservoir (although Coulomb blockade in 
its strict sense is absent). This leads to a suppression of 
the DOS, hence to a zero bias anomaly. Substituting Eq. 
(|l2|) in Eq. (||) one obtains for f^ 1 < r < 



S(t) S3 2r ln(0|r|) , 



(13) 



where 



R 



v<°) 



2nVL (2nh/e 2 ) 



(14) 



is the dimensionless resistance of the environment. This 



form of S leads to a power law DOS 
law cur rent- volt age characteristics [t 



and hence power 
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(El) <1 (to) <j 



FIG. |2|. Energy of the dot with N particles as function of 
the gate capacitance charge, q. Instantaneous tunneling of an 
electron into the dot, and subsequent slow relaxation of the 
charge are indicated by arrows, (a) Full screening, Co = oo. 
(b) Partial screening, a resonance scenario. 

Before concluding this section we draw a qualitative 
picture of the difference between the fully and partially 
screened interaction scenarios. Consider first the case 
of Cc — > oo and denote the number of electrons on the 
dot by TV. At long times the environment polarizes, such 
that the total neutralizing positive charge (including both 
the original positive background on the dot, Ao, and the 
accumulated charge on the gate capacitance, q) exactly 
balance off the charge N. Fig. ||a depicts the charg- 
ing energy of the dot as function of q, for various values 
of TV. This is given by E N {q) = e 2 {N - J\f - q) 2 /2C. 
The fact that the center of each parabola corresponds to 
En — is a manifestation of perfect screening. Upon 
adding an electron to the dot (N — > N + 1), the charging 
energy is initially increased, but at later times q relaxes 
towards the optimal value corresponding to N + 1. A 



where and Z^ are partition functions of uncou- 
pled dot and lead correspondingly, S(t) is defined by Eq. 
(D and 




The r.h.s. of Eq. ( |15[) may be considered as the 
grand-canonical partition function of a classical one- 
dimensional neutral gas consisting of positive and neg- 
ative charges which interact through the potential 
21nO|r| +S(t). The quantity 9, related to the resistance 
of the tunneling barrier by Eq. jl^ ) (cf. Eq. (Q)), rep- 
resents the fugacity of the gas. The interaction potential 
of the neutral gas deserves some elaboration. The first 
term, 21nfi|r|, came from the Green functions of non- 
interacting electrons in the dot and the lead. At finite 
temperature one should substitute r — > sui(itt/ (3)/3/tt. 
One should also understand that the interaction poten- 
tial given above is for large separation between charges 
(fir > 1), for t -y it should vanish. Note that S(t) 



very different situation arises for a partially screened in- 
teraction, Cq < oo. The charging energy is now given 
by E N (q) = e 2 (N - A/" - q) 2 /2C + e 2 q 2 /2C G . The cur- 
vature of each parabola corresponds to e 2 /2C + e 2 /2Ca, 
whereas the centers of the various parabola are located 
on the parabola e 2 /2(C + Cq)- Adding an electron to 
the system initially costs some energy (as in the previous 
case), but eventually the system relaxes to the partially 
screened state whose energy, although less than in the 
first moment after tunneling, is not in general equal to 
that of before tunneling event. By tuning Ao it is possible 
to have at most two states having the same final energy. 
Exactly this resonance scenario (with Ao = N + |) is 
depicted in Fig. ||b. This picture will be useful while 
considering multiple tunneling events, or quantum fluc- 
tuations of charge on the dot. 

III. QUANTUM FLUCTUATIONS OF CHARGE 

Our aim now is to address the full Hamiltonian, Eq. 
(^) , including the possibility of multiple tunneling events. 
To simplify calculations we first consider a partition func- 
tion of the system. We formally exactly expand it in 
powers of Ti. C oupUng and calculate separately each term 
of the expansion. Details of this procedure may be found 
in Appendix B. As a result one obtains the normalized 
partition function of the dot coupled to the lead in the 
following form 



(15) 



obeys this property as well as periodicity with (3 (by def- 
inition, cf. Eq. (|)). 

Before treating specific examples we make a few com- 
ments. Our gas which consists of the same number of pos- 
itive and negative charges is fundamentally different from 
the logarithmically interacting gas known from the the- 
ory of random matrices f2q| . The latter contains charges 
of the same sign placed in a confining potential. In our 
case the positions of the charges along the one dimen- 
sional (imaginary) time axis correspond to the times of 
the (instantaneous) tunneling events: positive charges - 
tunneling into the dot; negative - tunneling out of the 
dot. In general, any ordering of the charges, compatible 
with global neutrality is permissible. Neutrality follows 
from the fact that the partition function is given by a 
trace, thus forcing the initial and final charge states to 
coincide. At finite temperatures the interaction becomes 
periodic over the period < r < [3. The resulting gas is 
thus defined on a ring of circumference (3. 
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A. Fully Screened Interaction 

As has been discussed in Section E3, the fully screened 
interaction scenario corresponds to an infinite gate ca- 
pacitance, Cq — > oo, in the equivalent circuit (cf. Fig. 
|l|). This situation is equivalent to the case of a circuit 
with a single tunneling barrier, considered in Refs. 



in lowest order in tunneling. We provide here a system- 
atic treatment of higher order processes. Below we shall 
consider a purely Ohmic environment. The effective ac- 
tion to be substituted in the expression for the partition 
function is given (for (fi|rj — Tj\ » 1) by Eq. ([l]). Then 
the partition function is 



E 



q2M 



zmzw ^ M\M\ 

M=0 



f P ( 2M 

j. . .Id(n n ) . . . d(nr 2M ) exp I (2 + 2r) ^(-l)* 4 *' ln(r»|r 4 - Tj \ 

[ l< 3 



(17) 



This form breaks down when inter-charge distances 
are small, i.e. when the density of the Coulomb gas is 
high. This is a partition function of a classical Id gas of 
unit charges interacting through a 2d Coulomb (logarith- 
mic) potential at temperature (2 + 2r)~ 1 and fugacity 9. 
We denote this gas as non- alternating, meaning that any 
neutral sequence of charges is allowed. The correspon- 
dence between the quantum dot problem and the classi- 
cal Coulomb gas is summarized in a Table 1. Mapping 
onto the same problem has been found in the context of 
dissipative tunneling in a periodic potential |po|] and more 
recently in the study of the Luttinger liquid with an im- 
purity J2l| . The analogy stems from the observation that 
the Coulomb gas, Eq. (p~7|) , may be obtained from the 
partition function of a one dimensional quantum bosonic 
model with an action 



S[ip(x,r)] 



2r 



dxdT<p(d? + d?)<p 



-20Q 



/ dr cos (VA^ip{x = 0, t)) . (18) 



After integrating the bosonic field at all points except 
x = 0, one ends up with the effective action of the form 



S[<p(r)] = / dr 



1 



1 + r 



Lp\d T \tp — 29Q cos ( ViiTifi 



(19) 



Exponentiating the action and expanding in powers of 9 
one directly arrives at Eq. (17) where plays the role of 
the high energy cutoff. 

We follow here an approach based on a renormalization 
group (RG) study of the model (27j]2l]]. One integrates 
out the high energy (short time) degrees of freedom by 
continuously rescaling Q — > f2/£, where £ varies from 1 
up to fl/T ((3 = T^ 1 is the maximal possible size of the 
classical gas). As a result one arrives at the following 
well-known RG equations fl27|,pi| 



d9 
dr 



(20a) 
(20b) 



r 



The resulting RG trajectories on the (r, 9) plane are de- 
picted on Fig. |^a. In the context of a Luttinger liquid 
with an impurity, r represents the strength of the inter- 
action and is therefore a feature of the bulk system. The 
fugacity 9 represents the strength of the point-like im- 
purity. It is then clear that, while the interaction in the 
liquid may renormalize the impurity strength, the local 
impurity cannot renormalize the interaction strength in 
the bulk |^8). The separatrix r = corresponds to the 
presence of a Berezinskii-Kosterlitz-Thouless (BKT) like 
phase transition in a system of logarithmically interact- 
ing particles |2?|. The line r = corresponds to the 
non-interacting system (cf. Eq. (|l5|) with 5 = 0); evi- 
dently the tunneling transparency is not renormalized on 
this line. 




FIG. []. RG flows of the Coulomb gas: (a) a non-alternating 
gas; (b) an alternating gas. 

For positive r (which is the only physically meaningful 
scenario in the present context) the fugacity 9 ~ \W\ is 
renormalized down, 9 ~ £~ r . This provides a justifica- 
tion for a lowest order (in W) perturbation theory at low 
temperatures, which was employed in Section ||. Indeed, 
integrating Eq. ( |20a| ) and assuming that the relation be- 
tween fugacity and resistance, Eq. (|l6|), is conserved in 
the process of renormalization, one obtains 



(21) 



This zero-bias anomaly behavior has been derived previ- 
ously in a calculation which was non-perturbative 
in interaction, but perturbative in the tunneling. A sim- 
ilar power-law divergence of the resistance takes place in 
a repulsive Luttinger liquid with a single impurity pl| . 
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B. Partially Screened Interaction 



We now turn to a somewhat less trivial situation, that 
of a partially screened interaction. This corresponds to 
a finite Cq in Fig. [l| The ground state energies corre- 
sponding to different numbers of electrons in the dot are 
generically not degenerate. We refer though to a par- 
ticular scenario where the gate voltage (effectively Afo 
in Eq.(|l|)) is tuned such that the ground state energies 
corresponding to the charging states N, N + 1 are de- 
generate, cf. Fig. ||b. Injecting an electron to form 
any other charge state costs at least energy equal to 
il = e 2 /2(C + Cq) at long time scale. If we set the 
temperature to be less than this value one may neglect 
such states. Consequently the charge on the dot can 



fluctuate only between N and N + 1. The correspond- 
ing classical model is that of a gas of strictly alternating 
charges: a positive charge (tunneling into the dot) is fol- 
lowed by a negative charge (tunneling out of the dot). 
Similarly to the fully screened scenario discussed above, 
and following the same lines of derivation, the fact that 
the (two) available states are degenerate renders the long 
time (r ^> il^ 1 ) interaction logarithmic. The coefficient 
in front of the logarithm is determined by the slope of 
V(u>) at small u, which is given according to Eq. (noh by 



~] 2 



C + C c 



(22) 



As a result, the normalized partition function takes the 
form (see Appendix B) 



J 



2M / d(firi) / d{VlT 2 ) ... / d(QT 2M ) cxp I (2 + 2f) ^(-1) 4+J M^n - Tj\) 
M=0 n n n I *<J 



(23) 



The very same Coulomb gas was first discovered by 
Yuval and Anderson in the context of the Kondo prob- 
lem Q. In that case alternating positive and nega- 
tive charges correspond to the up and down spin flips. 
The mapping between an isolated dot near the Coulomb 
resonance and the anisotropic Kondo problem was first 
presented in Ref. 12 1. There the coupling to the en- 
vironment was absent (f = 0) and the inverse tem- 
perature of the gas was exactly 2 (see Table 1), ren- 
dering the corresponding Kondo problem to be always 
anti-ferromagnetic. Here we have introduced a coupling 
to the dissipative environment, which may change the 
situation qualitatively by driving the system through 
anti-ferromagnetic-ferromagnetic transition. Note that 
a small deviation from exact resonance conditions plays 
the role of Zeeman splitting due to a constant magnetic 
field in the Kondo model. In the case of the Kondo prob- 
lem both the fugacity and the inverse temperature are 
properties of the local spin, hence both may be renor- 
malized. The RG equations assume the form [ p^|j29| 

' W -r6 ; (24a) 



dlnf 

dr 
dlnf 



-40 2 



(24b) 



In the limit of extremely small fugacity, <C |f|, the gas 
is so dilute that the requirement of alternation does not 
play any role. Indeed in this case the RG equations of 
alternating gas coincide with Eqs. (20a) and (20b). To 



handle the case of finite fugacity one notices that 

r 2 - {29f = const . (25) 

As a result the RG trajectories are hyperbola in the (f, 9) 
plane. The RG diagram is depicted in Fig. ||b. The most 



striking feature of this diagram is the presence of the sep- 
aratrix at r — 29. In terms of the bare parameters of the 
problem the equation for the separatrix has the form 



TT^/nr = 1 . 



(26) 



The low temperature behavior of the system is qualita- 
tively different, depending on whether the bare parame- 
ters are chosen to be below or above the separatrix line. 

Let us first consider the system below the separatrix, 
rt < (irf)~ 2 . In this case the fugacity scales down, and 
the resulting Coulomb gas is dilute (the system is in 
the "cold" phase, i.e., overdamped by the environmen- 
tal resistance). Lowest order perturbation theory should 
be adequate while renormalization provides some correc- 
tions. One can easily integrate Eqs. ( |24a| ) and (24b) by 
putting 

f(0 = Acoth(Aln£ + B); 20(f) = Aj sinbLAlnf + B) , 

(27) 



where A = 
one obtains 



=2 



(20) 2 and cothB = f/A. For T < f2 



r*(T)~T 



-2 A 



(28) 



If <C f , then A k f and Eq. ( |28| ) coincides with the 
result of perturbation theory (in the tunneling strength), 
cf. also Eq. (|2l]). In the more general case the low tem- 
perature exponent is somewhat smaller than the naive 
perturbation theory prediction, indeed A < r. 

At the separatrix line, Eq. (|26|), the system undergoes 
a zero-temperature BKT phase transition. Above the 
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separatrix line, r t > (7rr)~ 2 , the fugacity exhibits reen- 
trance behavior and eventually is renormalized up. The 
corresponding Coulomb gas is in the "hot" or plasma 
phase. Note that the dot without a dissipative environ- 
ment (f = 0) is always in the hot phase. This latter case 
has been extensively studied since the work of Matveev 
[ p2[ (see e.g. Ref. Jl3|]). The temperature below which 
the reentrance behavior takes place is analogous to the 
Kondo temperature, Tk- To find Tk we integrate Eqs. 
(|24a|) and (24b) by the substitution 



f(f) = -itan(iln^-B); 20(£) = A/cos(Aln£ - B) , 

(29) 

where A 2 = (29) 2 —r 2 and tanf? = r/A. The Kondo tem- 
perature may be estimated by equating the argument of 
the tan to tt/2. This way one obtains 



Tk = fl exp < — 



(30) 



The Kondo temperature is maximal for an isolated dot, 
f = 0. In this case Eq. © reduces to @ T K = 
exp{— 7r 2 y / ?t'/2} . Thus even in the best case the 
Kondo temperature is exponentially small in the root of 
the tunneling resistance. In order to assure the validity 
of this theory the inequality Tk > A should be satisfied, 
otherwise the assumption about a continuous spectrum 
("metallic dot") is not valid and one has to deal with the 
discrete spectrum limit. For a finite value of the envi- 
ronmental resistance, f, the Kondo temperature rapidly 
approaches zero once the parameters of the system come 
closer to the separatrix line, Eq. (|26|). This makes the ex- 
perimental observation of the reentrance behavior quite 
difficult in practice. In contrast, the "cold" phase behav- 
ior described above seems to be relatively easy to observe 
by measuring the temperature dependence of the heights 
of the Coulomb blockade peaks. 



IV. SUMMARY AND DISCUSSION 

We have addressed here the physics of a (quasi)zero- 
dimensional quantum dot weakly coupled to an ideal lead 
and capacitively (and not necessarily weakly) coupled 
to an external dissipative environment. The fact that 
charges on the dot interact with the environment which 
has its own dynamics is translated into an effective time 
dependent potential felt on the dot. In other words, when 
a charge is added to the dot it polarizes the environment. 
This process is not instantaneous — it is characterized 
by an effective RC time. The delayed response of the 
environment leads to a non-trivial effective time depen- 
dent interaction term. The (weak) coupling of the lead 
to the dot means that charges can tunnel in and out of 
the dot, resulting in quantum fluctuations of the charge 



on the dot. We have considered a single channel cou- 
pling. For this to be the case we require that the linear 
dimension of the opening connecting the dot to the lead 
does not exceed the Fermi wave length. The methods we 
employed (mapping onto the Coulomb gas and pertur- 
bative RG) are restricted to the case of relatively weak 
tunneling coupling, that is barrier resistance larger than 
2irh/e 2 . Although we did not restrict ourselves to any 
finite order in the tunneling coupling, a consistent treat- 
ment of the strong coupling fixed point is outside the 
scope of this article. The nature of an appropriate fixed 
point was investigated in Ref. ]30| ]. Our main message 
here is rather the existence of a separatrix in the param- 
eter (resistance) space of the problem, which separates 
the basins of attraction of the weak and the strong cou- 
pled fixed points. 



g 

1/T 



FIG. y. Schematic behavior of the conductance (in units 
e 2 /h) as function of inverse temperature (in units 1/fi) for 
the overdamped (full line) and underdamped (dashed line) 
scenario; r = .75. 

We have proceeded by mapping our quantum model 
onto a classical one-dimensional neutral Coulomb gas. 
The nature of this classical gas depends crucially on 
whether the Coulomb term is fully or partially screened. 
In either case we obtain the corresponding renormaliza- 
tion group flows. In the former case (corresponding to 
the single tunnel barrier connected to a dissipative cir- 
cuit J|,[llJ ) the flow is towards the "cold" phase (weak 
coupling fixed line; Fig. ||a). The results of the pertur- 
bative treatment are naturally reproduced in this 

case. That is, at low temperatures the tunneling coupling 
is suppressed as a certain (positive) anomalous power of 
temperature (cf. Eq. (Bj)). The case of the dot with a fi- 
nite residual capacitance (partial screening) is more intri- 
cate. In particular, we found, in agreement with Matveev 
]l2[ and others (l3) that for sufficiently small dissipation 
of the environment the problem is equivalent to a anti- 
ferromagnetic Kondo. The RG flow is towards the strong 
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coupling, characterizing by the ideal transparency of the 
contact. If, however, the resistance of the environment 



exceeds a certain critical value (cf. Eq. (26)) determined 
by the bare tunneling resistance, the situation changes 
dramatically. The system turns to be overdamped and 
its low temperature behavior is determined by a weak 
coupling fixed line, similar to those of the fully screened 
situation (with a slightly modified exponent). 

The physics of this damping of charge fluctuation by 
the environment is simple. Consider a dot which is tuned 
by the gate electrode into a resonance, i.e., the cost (in 
free energy) of removing one electron from the leads and 
adding it to the dot is nearly zero. (We presently assume 
no voltage difference between the two leads connected to 
the dot). The environmental resistance produces noisy 
voltage on the gate electrode. Thus, at any given mo- 
ment, the dot is slightly driven out of resonance. Not 
surprisingly, charge fluctuations are suppressed by this 
(time-dependent) deviation from the exact resonance 
condition. What is more surprising is an existence of a 
sharp transition at zero temperature (in fact at T < Tk 
) between two regimes as a function of the environmen- 
tal resistance. At a finite temperature one has a smooth 
crossover between these two limits. The schematic tem- 
perature behavior of the tunneling conductance for the 
two cases distinguished by the value of an environmental 
dissipation is depicted on Fig. |. To be within a scenario 



of an alternating Coulomb gas the initial temperature 
should be less than both e 2 /2(C + C G ) and (RCy 1 . 

Throughout the above analysis it has been assumed 
that the relation between the resistance and the fugacity, 
Eq. ( |l6| ) , remains valid under the renormalization group 
transformation. It is possible, though, to establish a di- 
rect relation between the Coulomb gas description and 
the conductance. We first note that the current oper- 
ator, J, is given by the time derivative of the number 
operator N = a i a a 



k a 

The linear conductance may be expressed, employing the 
Kubo formula, through current-current correlation func- 
tion. Evidently, the latter is expressible through the two 
following correlation functions 



r 2 ) 

T2) 



(TW(n)W + (T 2 )) , 

(TW{ n )W{T 2 )) , 



(32a) 
(32b) 



where T denotes an (imaginary) time ordering. Both F 
and F may be naturally interpreted in a Coulomb gas 
language. Following the same procedure depicted in Ap- 
pendix B we arrive at the following expression for e.g. 

F{ n - r 2 ) 



J 



1 oo 

F( n -r 2 ) = -J2 



Z ^ M\M\ 

M=0 



^ ^ ( 2M+2 

/... d((lT 3 )---d(nT 2M+2 )expt (-l) l+J [21n(r!|T i -r,|) + 5'(r. 1 
o I *<i 



(33) 



r 



This is the partition function of a one-dimensional 
classical neutral gas, the only difference is that now two 
charges of opposite charge are frozen at the points t\ 
and T2 respectively, all others (whose number is not 
fixed) are free to move. We may write F{t\ — t 2 ) oc 
exp{— U{\t\ — t 2 \)}, where U is the effective screened in- 
teraction between the two frozen opposite charges. In an 
analogous way F may be expressed through the screened 
interaction of two positive (negative) charges. It is im- 
portant to stress that a single particle propagator (more 
generally: a correlator of an odd number of particles) 
does not have a gas equivalent. The reason is that dif- 
ferent configurations of "charge" correspond to terms of 
different signs (an exchange of two charges gives rise to 
a sign change). Thus not all terms can be assigned a 
positive weight. 

There are various ways to generalize and extend the 
model discussed here. First, we may include interactions 
in the lead and even between electrons in the lead and 
electrons on the dot. In an obvious notation we denote 
the two types of interactions by Vu and Vqi respectively, 
while the interaction associated with the dot is Vqo- We 



assume that these (effective) interactions may be time de- 
pendent, but we do not include space dependence (only 
capacitive component). The only way these interactions 
modify the results is through the definition of the effec- 
tive action, Eq. (0), where one substitutes 



V{lo) -» V Q0 + V H - 2V i . 



(34) 



To see why this is the case one should proceed exactly 
along the lines of Appendices A and B, introducing a 
doublet of Hubbard-Stratonovich fields, 4>o-,4>i- 

Inclusion of the spin degree of freedom is trivial (re- 
definition of the DOS) as long as the tunneling does 
not conserve spin. If spin is conserved in the process 
of tunneling this is translated into a gas with charges 
(positive and negative) of two "colors" . Charges of the 
same color interact among themselves via the potential 
2 In fl | t | + S (r ) , as in Eq. (|l^) . Charges of different colors 
interact through the potential S(t). This implies that in 
the absence of real interaction in the dot (S = 0), the 
gas factorizes into two mutually non-interacting gases. 
Using an analogy with the multichannel Kondo problem 
one may claim that the separatrix line still exists in the 
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system's phase diagram. Its expression through the bare 
parameters of the dot requires, however, some more cal- 
culations. 
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APPENDIX A: SINGLE PARTICLE GREEN FUNCTION 



Let us hrst calculate the imaginary time single particle Green function of an isolated dot described by a Hamiltonian 
Eq. (|l|). It may be written as |u| 



Gafa - Tf,fl) 



Z{n) 



V[a a (T)a a (T)]e 



-S[a a ,a a ] 



a a (n)a a (Tf), 



with the fermionic action given by 

P 



S[a a ,a a ] = J dr 



ut°] 



^ a a (t) (d T +e a ~ n)a a fa) + —— 



^2a a (r)a a (T) - JV 



(Al) 



(A2) 



here Z(n) is the partition function and /x is the chemical potential. Splitting the interaction term in the action by 
means of the Hubbard-Stratonovich transformation with the auxiliary Bose held, c/>(t), one obtains 



(r)]e ° 



- jf dT[0(r)[2y[°]]- 1 0(r)- l A^0(r)] 



— J dr a a (d T +e a — f.L+i4>{r))a a 

V[a a a a ]e a a fai)a a faf) 



r -/dr[0(T)[2y[°)]- 1 0(T)-iM^(T)] 



Z(jj) 



(A3) 



with the same transformations in Z(/i). Here Z^(n) and Go" fat, r/, /i) are respectively the partition and Green 
functions of non-interacting electrons in the time dependent (but spatially uniform) potential, i4>(r). These quantities 
may be easily calculated, taking advantage of a gauge transformation. The result is 



ZM(/i) = Z[°](/i-^ ) 



(A4) 



(A5) 



where Z^ \/j.) = exp{— (//)} and Ga\e n , /i) = (ie n — e a + fi) 1 are the partition and Green functions of non- 
interacting electron gas; we have introduced Matsubara representation for the boson field, 4>{t): 



<t>m = 1 / (1t4>(t) exp{iw m r}, 
Jo 



uj m = 2irmT. For latter convenience we rewrite the exponent in Eq. flAq ) in the following form 

expli dr[0(r) - 4> Q ] > = exp I (3 ^ 



jTi\Tf 

-m'->ra 



(A6) 
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where Jm' r/ is the Matsubara transform of the following function 

J? lTf =5(T-Ti)-S(T-T f ). 

Transforming the functional integral over </>(r) to integrals over the Matsubara components, 4> m , we obtain 



(A7) 



Qa(n-Tf) = 



1 



ZQi) 



b e 



fopV^-^o-iclioAfa+n^Hp-iMgiO]^. _ Tf)fl - i(j) Q ) 



/ Yl d( f' m ex p \ p H 

m^O ! m^O 



jTi\Tf 



2V^ 



(A8) 



For large enough systems (A -C T) the integral over the static component, tfio, may be calculated in a saddle 
point approximation leading to Q^iji — Tt,~p), where the stationary point, /I, is the real solution of the equation 
(fx — ~p) /V^ + N a + dfl^ (~p) /dfi = 0. The remaining integrals (over <p m for m 7^ 0) are purely Gaussian. As a result 
one obtains 



where 



m#0 m 



(A9) 



(A10) 



We shall include now the environment, Eq (|^), in the calculations of the Green function. Assuming that the noise 
is Gaussian with zero mean value, one first perform the averaging over noise realizations as 



(...), 



-3 // ' dTdT , r)(T)K- 1 (T-T')r](T') 

V[r](r)]e ° 



(All) 



where K(t — t') = (^(tWt')). An important observation is that the partition function of the dot, Z(fi), is not 
affected by the noise term |3^] (cf. Eq. (A4)). Averaging Eq. ( |Al| ) over the Gaussian noise, Eq. ( All ), leads to an 
effective fermionic action with interaction which is non-local in time, 



Sint [a a ,a a ] =2 drdT ' ® aaa (^ [0l ^( r _ r ') _ K ( T ~ r ')) ® aCLa - 



(A12) 



As a result, one obtains an effective renormalization (screening) of the zero-mode interaction potential 

VM - V{uj m ) = VM - K(uj m ). (A13) 
Further calculations follow exactly the same steps as outlined above. The final result is given by Eqs. (0) and (@). 



APPENDIX B: MAPPING ONTO A COULOMB GAS 

Our aim here is to derive an expression for the partition function corresponding to the full Hamiltonian, Eq. (|(j). 
Performing averaging over the noise and the Hubbard-Stratonovich transformation as described in Appendix A one 
obtains 



2?M«p - \Z[<t>] (Bl) 



where 
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Z[4>] = I V[aadd] exp I - [ drC \ , (B2) 



C = ^2a a (d T + e a +i<f>(r))a a + d k (d T + e k )d k + W^d k a a + W*^a a d k . (B3) 



ka kc 



One may integrate out now all fcrmions of the dot and lead except the two living at the point of the tunneling contact 
by introducing the following resolutions of unity 

1 = J V[s]S(s - O = / °] ex P I i J dr °( s -J2 aa U ' ( B4 ) 
a " I o J 

and similarly 

1 = J V[r, p] exp |i J drp(r - £ d k ) 1 , (B5) 

where s, r, a, p are Grassman variables. We also introduce the two corresponding identities for the conjugate variables. 
Performing the Grassman integration over the original variables {a, a, d, d} and then over the auxiliary {a, a, p, p} one 
obtains 

Z[<f>] = Z [a] Z [l] (^expi^- JdT(Wfs + W*sr)^J , (B6) 
where the angular brackets denote integration over fields s and r with the following (^-dependent measure: 

det 0^ dot Jv[ssfr] exp i~f f drdT ' + f^ 1 ' 1 )" 1 ?-] 1 . (B7) 

Here we have introduced notations for the traces of the Green functions of the dot and the lead respectively 

exp{/?£ 



g® (r, r>) = £ QM (t, r>) = „M \^_mj^_ 1 

L — ' T — T 

a 

gW iT -T>) = ^gV(T-T>)=vW_L_. (B9) 

k 

We have assumed that both the dot and the lead have continuous spectra, characterized by constant density of states 
and i/M respectively), and employed Eqs. (|A5|)-(A7). At finite temperatures one has to substitute 



r — t' sin7r(r — r')//3 



11/(3 (BIO) 



The next step is to expand the exponent in Eq. (B6) to infinite order in W, yielding 
Z[<f>] _ ^ \W\ 2M (2M)I 



p p^-pejw (2M)\ f f 

Z [0] Z H) " £ (2M)! Ml Ml J ■••y dr i--- dT 2Af(sis 2 ... S2M _ 1 s 2M )(f 1 r 2 ...f 2M _ 1 r 2M ) , (Bll) 

M=0 

where Si = s(t,), i = 1,2 . . . 2M. The combinational factor (2M)!/(M!M!) represents the number of possibilities to 



rearrange time sequences in the angular brackets. Due to the Gaussian nature of the measure, Eq (B7), one may 
apply Wick's theorem, resulting in the M x M Slatter determinant construction 
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(nr 2 ■ ■■r 2 M-ir 2 M) = det |<? M (t 2 „_i - r 2m )| 



(B12) 



n, m — 1, 2 . . . M. We employ now the explicit form of the Green function, Eq. (B9), along with the known properties 
of a Cauchy determinant |l7]] to write 



II ( T 2n — T 2m ) II ( T 2n-1 — T 2m _l) 

{r\r 2 . . .r 2M -ir 2M ) = -{v u ) TT1 \ 

11 \ T 2n - T 2m -\) 



(B13) 



In the same manner 



(siS 2 . . . S2M-1S2M) = ~(v^) M — 



Y[ ( T 2,i-T 2m ) Y\ ( r 2, 1 -l-T 2m -l) 



{T2n-T 2m -l) 



ri( 

2--- T 2M \ T 1 '- T 2M-1 



where J m is the Fourier transform of the imaginary time function 

2M 

jr 2 ...r 2M |r 1 ...T 2 M-l = S{t - Tj 

3 = 1 



(B14) 



(B15) 



The exponent in Eq. ( B14 ) is the only part of Z[<fi] which depends on the Hubbard-Stratonovich field , 4>{t). Thus 
integration over the field (/>, Eq. (Bl), may be easily performed for each term in the sum (Eq. (Bll)) separately. 
Indeed 



|exp 




■'—m y vm 



i.Jrr 



2V(uo m ) W77 



= exp 




(B16) 



One may now return to time representation and employ the definition of J T , Eq. ( B15 ), as well as Eq. (|8|) 

2M 

to write the exponent on the r.h.s. of Eq. ( B16| ) as (— l) !+J S'(T i — Tj). Note also that the product of the 

i<j 

two Cauchy determinants (Eqs. ( B13| ) and ( B14| )), being positive definite, may be written as an exponent of 

2M 

Y^, ( — 1) I+J 2 In \t.i — Tj\. Finally one obtains Eq. ( |l5| ) for the partition function of the dot connected to the lead. 

i<j 

I 
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